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We study the real-time domain-wall dynamics near a quantum critical point of the one-dimensional 
anisotropic ferromagnetic spin 1/2 chain. By numerical simulation, we find the domain wall is 
dynamically stable in the Heisenberg-Ising model. Near the quantum critical point, the width of 
the domain wall diverges as (A — 1) _1//2 . 
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I. INTRODUCTION 



Recent progress in synthesizing materials that contain 
ferromagnetic chains^^ provides new opportunities to 
study the quantum dynamics of atomic-size domain walls 
(DW). On the atomic level, a DW is a structure that is 
stable with respect to (quantum) fluctuations, separating 
two regions with opposite magnetization. Such a struc- 
ture was observed in the one-dimensional C0CI2 • 2H 2 
chain£i&. 

In an earlier paper—, we studied the propagation of spin 
waves in ferromagnetic quantum spin chains that sup- 
port DWs. We demonstrated that DWs are very stable 
against perturbations, and that the longitudinal compo- 
nent of the spin wave speeds up when it passes through a 
DW while the transverse component is almost completely 
reflected. 

In this paper, we focus on the dynamic stability of the 
DW in the Heisenberg-Ising ferromagnetic chain. It is 
known that the ground state of this model in the subspace 
of total magnetization zero supports DW structures^. 
However, if we let the system evolve in time from an 
initial state with a DW structure and this initial state is 
not an eigenstate, it must contain some excited states. 
Therefore, the question whether the DW structure will 
survive in the stationary (long-time) regime is nontrivial. 

The question how the DW structure dynamically sur- 
vives in the stationary (long-time) region is an interesting 
problem. In particular, we focus on the stability of the 
DW with respect to the dynamical (quantum) fluctua- 
tions as we approach the quantum critical point (from 
Heisenberg-Ising like to Heisenberg). We show that the 
critical quantum dynamics of DWs can be described well 
in terms of conventional power laws. The behavior of 
quantum systems at or near a quantum critical point is 
of contemporary interest^. We also show that the DW 
profiles rapidly become very stable as we move away from 
the quantum critical point. 



II. MODEL 



The Hamiltonian of the system is given by 
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N-l 

H =- J H ( SX n S n+l + S l S l + l + AS*S* + 1 ), (1) 
n=l 

where N indicates the total number of spins in the spin 
chain, and the exchange integrals J and J A determine 
the strength of the interaction between the x, y and z 
components of spin 1/2 operators S„ = S% , S*). 
Here we only consider the system with the ferromagnetic 
(J > 0) nearest exchange interaction. It is well known 
that |A| = 1 is a quantum critical point of the Hamil- 
tonian in Eq. |T]), that is, the analytical expressions of 
the ground state energy for 1 < A and — 1 < A < 1 are 
different and singular at the points A = ±1—. 

In Refj££ Gochev constructed a stable state with DW 
structure in both the classical and quantum treatments of 
the Hamiltonian (fTJ) . In the classical treatment, Gochev 
replaces the spin operators in Eq. fT]) by classical vectors 
of length s 

S* = s cos 6 n , S„ = s sin 6 n cos ip n , S% = s sin n sin ip n , 

and then uses the conditions SE/S9 = and ip n = const. 
to find the ground state. In the ground state, the mag- 
netization per site is given by^ 



= stanh(n — no)a, 
S% = s cos if sech(n — no)a, 
S% = ssincp sech(n — no)a, 



where 



a = ln[A + VA 2 -!], 



(2) 



(3) 



ip is an arbitrary constant, and no is a constant fixing the 
position of the DW. The corresponding energy is 



E DW = 2s JA tanner. 



(4) 
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In the quantum mechanical treatment, Gochev first con- 
structs the eigenfunction of a bound state of k magnons 9 

\tpk) = A n ^■>n 1 m 2 ...m k S mi S Jn2 ...S Jnk |0) , (5) 

{mi} 



where 



B 



= cosh(z — l)u/ cosh(za), 



(G) 



(7) 



nnnm^iUiiiui 

Mill II 11 + UUUU 

(b) 

MIIIIII'^UUUU 



(c) 



nnnninjjjjiiui 

lllllllii'niiuuu 

ffiffiffiflMHIil 



A- 2 =n^/(i-^), 



and the corresponding energy is given by2. 



1 



(~k = — JAtanhcrtanh ka. 



(8) 



(9) 



Then he demonstrated that for the infinite chain, the 
linear superposition 



oo , 1 

l^no) = A CXP 1 ~0 a 



I 2 a 



(10) 



where 



n = N + a, \a\ < 1/2, N -> oo, (11) 



A- 2 = 



'j — — OO 
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(12) 



is the quantum analog of the classical domain wall, in 
which (S*) , (S*) , (S%) are given in Eq. ©, and the en- 
ergy coincides with Eq. (fj|. 

Gochev's work confirmed the existence of the DW 
structures in the one-dimensional ferromagnetic quantum 
spin 1/2 chain. In the infinite chain, the exact quantum 
analog of classical DW is represented by |0„ o ). In the fi- 
nite chain, the DW structure exists as a bound fc-magnon 
state \ipk)- The main difference between these two states 
is the distribution of magnetization in the XY plane. In 
the infinite chain, the change of the magnetization oc- 
curs in three dimensions, according to Eq. @, but in 
the finite chain (S%) — (S%) — for all spins. 

Now we consider (S*) of the bound state iV'fc) in the 
case that the number of flipped spin is half of the total 
spins, i.e., k = N/2 and N is an even number. Even 
though the formal expression for \ipk) is known, the ex- 
pression for (S!//) in this state (for finite and infinite 
chains) is not known. For finite N, the ground state in 
the subspace of total magnetization M = can, in princi- 
ple, be calculated from Eq. ((?}]). However, this requires a 
numerical procedure and we loose the attractive features 



FIG. 1: The magnetization {S^} in the ground state of the 
subspace of total magnetization M = 0, generated by the 
power method. The parameters are: (a) A = 1.05, (b) A = 
1.1, (c) A = 1.2, (d) A = 2. The total number of spins in the 
spin chain is N = 20. It is clear that there is a DW at the 
centre of the spin chain. Furthermore there is no structure in 
the XY plane, that is, (S*) = (S*) = 0. 
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FIG. 2: Top picture (a): Initial spin configuration at time 
t/r — 0; Bottom pictures (b,c,d,e,f,g,h,i): Spin configuration 
at time t/r = 100; Bottom left pictures (b,c,d,e): DW struc- 
tures disappear or are not stable. The parameters are: (b) 
A = (XY model), (c) A = 0.5 (Heisenberg-XY model), (d) 
A = 1 (Heisenberg model), (e) A = 1.05 (Heisenberg-Ising 
model); Bottom right pictures (f,g,h,i): DW structures are 
dynamically stable in the Heisenberg-Ising model. The pa- 
rameters are: (f) A = 1.1, (g) A = 1.2, (h) A = 2, (i) A = 20. 
The total number of spins in the spin chain is N = 20. 



of the analytical approach. Indeed, it is more efficient to 
use a numerical method and compute directly the ground 
state in the subspace of total magnetization M = 0. In 
Fig. HI we show some representative results as obtained 
by the power method^ for a chain of N = 20 spins. In 
all cases, the domain wall is well-defined. Obviously, be- 
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TABLE I: The energy E = JA/2 of the initial state |$> (see 
Fig. [2ja)) and the ground state Eg in the M — subspace, 
both relative to the ground state energy of the ferromagnet. 



A 


E 




£(18) 


p(20) 


F (22) 


1.05 


0.53 


0.16 


0.16 


0.16 


0.16 


1.1 


0.55 


0.23 


0.23 


0.23 


0.23 


2 


1.00 


0.87 


0.87 


0.87 


0.87 


5 


2.50 


2.45 


2.45 


2.45 


2.45 



cause we are considering the system in the ground state, 
the magnetization profile will not change during the time 
evolution. 

To inject a DW in the spin chain, we take the state |$) 
with the left half of the spins up and the other half down 
as the initial state (see Fig. [2ja) for TV = 20). The state 
|$) corresponds to the state with the largest weight in the 
bound state \ipk) with k = N/2, because \B mim2 ,„ mk | 2 
reaches the maximum if m, = i for all i = 1,2, ..,N/2 
(note \vj,\ < 1). It is clear that |$) is not an eigenstate of 
the Hamiltonian in Eq. ([T]). The energy of |$) , relative to 
the ferromagnetic ground state, is JA/2, and its spread 
(($\H 2 \<f>)-($\H\$) 2 y/ 2 = J/2. In Tablei we list some 
representative values of the energy in the initial state (see 
Fig. [DJa)) and in the ground state of subspace M = 
(see Fig. Q}. 

A priori, there is no reason why the DW of the ini- 
tial state |$) should relax to a DW profile that is dy- 
namically stable. For A ~ 1, the difference between en- 
ergy of the initial state and the ground state energies 
for N = 16, 18, 20, 22 is relatively large and the relative 
spread in energy (1/A) is large also, suggesting that near 
the quantum critical point, the initial state may contain 
a significant amount of excited states. Therefore, it is 
not evident that a DW structure will survive in the long- 
time regime. In fact, from the numbers in Table I, one 
cannot predict whether or not the DW will be stable. 
For instance, for A = 1.05 and N = 16, 18,20, the DW 
is not dynamically stable whereas for N = 22 it is stable 
but the energies (see first line in Table I) do not give a 
clue as to why this should be the case. On the other 
hand, by solving the time dependent Schrodinger equa- 
tion (TDSE), it is easy to see if the DW is dynamically 
stable or not. 



III. DYNAMICALLY STABLE DOMAIN WALLS 

We solve the TDSE of the whole system with the 
Hamiltonian in Eq. |T]) and study the time-evolution 
of the magnetization at each lattice site. The numeri- 
cal solution of the TDSE is performed by the Chebyshev 
polynomial algorithm, which is known to yield extremely 
accurate independent of the time step use d 15 i 16 i 17 i 18 . We 



adopt open boundary conditions, not periodic bound- 
ary conditions, because the periodic boundary condition 
would introduce two DWs in the initial state. In this pa- 
per, we display the results at time intervals of r = 7r/5 J, 
and use units such that K = 1 and J = 1. 

The initial state of the system is shown in Fig. [2fa). 
The spins in the left part (n = 1 to 10) of the spin chain 
are all "spin- up" and the rest (n = 11 to 20) are all 
" spin-down" . Here " spin- up" or " spin-down" correspond 
to the eigenstates of the single spin 1/2 operator S£. 

Whether the DW at the centre of the spin chain is 
stable or unstable depends on the value of A. In Fig. 
[5{b,c,d,e,f,g,h, and i), we show the states of the system 
as obtained by letting the system evolve over a fairly long 
time (t = 500J/7r). It is clear that the DW totally disap- 
pears for < A < 1, that is, in the XY, Heisenberg-XY 
and Heisenberg spin 1/2 chain, the DW structures are 
not stable. For the Heisenberg-Ising model (A > 1), the 
DW remains stable when t > 500 J /ir (see Ref.— ), and its 
structure is more sharp and clear if A is larger, so we 
will concentrate on the cases A > 1. One may note that 
the values of A in Fig. [2je,f,g,h) are the same as in Fig. 
[T][a,b,c,d), but that the distributions of the magnetiza- 
tion are similar but not the same. This is because the 
energy is conserved during the time evolution and the 
system, which starts from the initial state shown in Fig. 
[Ha), will never relax to the ground state of the subspace 
with the total magnetization M = 0. 

In order to get a quantitative expression of the width 
of DW, we first introduce the quantity S* (ti, taj A) (n = 
1,2,..., N) as the time average of the expectation value 
(S£ (t)) of nth spin: 



S*(ti,t 2 ;A) 



t 2 -h 



(13) 



We take the average in Eq. (|T3|) over a long period dur- 
ing which the DW is dynamically s table. In Fig. [3J we 
show some results of S*(ti,t 2 ;A) for the Heisenberg- 
Ising model, where we take t\ — 101r, t 2 — 200t and 
various A. We find that each curve in Fig. [3] is symmet- 
ric about the line n = (N + 1) /2, and can be fitted well 
by the function 



S% {t\,t 2 ; A) = oa tanh 



n-(JV + l)/2 



(14) 



The values of A we used and the corresponding values 
of a A> are shown in Table HH As we mentioned 
earlier, Gochev^ constructed an eigenstate of the one- 
dimensional anisotropic ferromagnetic spin 1/2 chain in 
which the mean values S*, and <S^ coincide with the 
stable DW structure in the classical spin chain, that is 



1 



(15) 



where no is the position of the DW (in our notation, this 
is (N + 1) /2). The fitted form of S* {h,t 2 ;A) in Eq. 
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FIG. 3: (Color online) S*(ti,t2',A) as a function of n for 
different A. Here ti = lOlr, ti = 200r. We show the data 
for A = 1.05, 1.06, 1.1, 1.2, 1.3, 1.5, 2, 5 and 20 only. The total 
number of spins in the spin chain is N = 20. 



(fT¥jl is similar to Eq. (Ti~5"]) . From Table HT1 it is clear that 
as A increases, |oa| converges to 1/2, in agreement with 
Eq. (|T5|) . From the comparison of 6a and 1/a in Fig. 
21 it is clear that the dependence on A is qualitatively 
similar but not the same. This is due to the fact that 
Gochev's solution is for a DW in the ground state whereas 
we obtain the DW by relaxation of the state shown in Fig. 
Ha). 

We want to emphasize that the meaning of 
S* (h,t 2 ;A) in Eq. ^ is different from {S*} in Eq. 
(fT5| . The former describes the mean value of (S* (t)) 
in a state with dynamical fluctuations, while the latter 
describes the distribution of (S*) in an exact eigenstate 
without dynamical fluctuations. 

Next we introduce a definition of the DW width. From 
Eq. (|14[) , we can find n\ and n 2 which satisfy 



S*, (ti,t2;A) = l/4, 
-1/4, 



S z n2 (h,h;A) 



(16) 



that is, when S% (t\, t 2 ] A) equals half of its maximum 

value (1/2). Here ri\ and n 2 are not necessarily integer 
numbers. Now we can define the DW width W as the 
distance between ri\ and n 2 : 



W = \ni - n 2 \ 



(17) 



Clearly, the width of the DW becomes ill-defined if it 
approaches the size of the chain. On the other hand, the 
computational resources (mainly memory), required to 
solve the TDSE, grow exponentially with the number of 
spins in the chain. These two factors severely limit the 
minimum difference between A and the quantum critical 
point (A = 1) that yields meaningful results for the width 
of the DW. Indeed, for fixed N, A has to be larger than 



TABLE II: The values of A we used in our simulations and 
the corresponding a a, &a fitted by Eq. (|14p for a spin chain 
of N = 20 spins. 



A 




b A 


A 


OA 


b A 


1.05 


— 0.263 


3.659 


1.8 


— 0.493 


0.524 


1.06 


— 0.330 


3.171 


1.9 


— 0.494 


0.488 


1.07 


-0.377 


2.850 


2 


-0.495 


0.460 


1.08 


-0.406 


2.673 


2.1 


-0.495 


0.436 


1.09 


-0.424 


2.534 


2.2 


-0.496 


0.416 


1.1 


-0.435 


2.396 


2.5 


-0.497 


0.370 


1.15 


-0.462 


1.996 


3 


-0.498 


0.322 


1.2 


-0.471 


1.626 


4 


-0.499 


0.270 


1.25 


-0.476 


1.330 


5 


-0.499 


0.240 


1.3 


-0.479 


1.142 


6 


-0.500 


0.220 


1.35 


-0.481 


0.959 


7 


-0.500 


0.206 


1.4 


-0.483 


0.869 


8 


-0.500 


0.195 


1.45 


-0.485 


0.770 


9 


-0.500 


0.187 


1.5 


-0.487 


0.719 


10 


-0.500 


0.179 


1.6 


-0.489 


0.629 


15 


-0.500 


0.156 


1.7 


-0.491 


0.568 


20 


-0.500 


0.141 
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FIG. 4: Comparison of 6a and 1/<t as a function of A. The 
total number of spins in the spin chain is N = 20. 



the "effective" critical value for the finite system in order 
for the DW width to be smaller than the system size. 
Although the system sizes that are amenable to numeri- 
cal simulation are rather small for present-day "classical 
statistical mechanics" standards, it is nevertheless possi- 
ble to extract from these simulations useful information 
about the quantum critical behavior of the dynamically 
stable DW. 

In Fig. O we plot W as a function of A (1.06 < A < 
20). By trial and error, we find that all the data can be 
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As a check on the fitting procedure, we apply it to 
the data obtained by solving for the ground state in the 
M = subspace. In view of Eq. (fl~3|) and (fT4|) . we may 
expect that Eq. (fl~8"|) fits the data very well and, as shown 
in Fig. [8j this is indeed the case. 

If we fit the data to 



W (A) = Wo (A - A 



N) 



(20) 



without assuming a priori value C, we find that C de- 
pends on the range of A that was used in the fit, as 
shown in Fig. [9] Remarkably, we find that C w 0.57 if 
we fit the data for a large range of A's and that C ap- 
proaches 1/2 if we restrict the value of A to the vicinity 
of the critical point. 



FIG. 5: The DW width as a function of A in a spin 
chain of TV = 20 spins. The black dots are the sim- 
ulation data and the solid line is given by W (A) = 

A N /\n^A - e N ± [(A - e N f - l] 1/2 | + B N with e N = 

0.046 ± 0.001, A N = 2.16 ± 0.06 and B N = -0.485 ± 0.068. 



TABLE III: The values of e N , A N and B N in Eq. JT8J) for a 
spin chain of N = 16, 18, 20, 22 and 24 spins. For the fits, 
we used all the data for A < 5. 



N 




A N 


Bn 


16 


0.065 ±0.001 


2.08 ±0.10 


-0.493 ±0.142 


18 


0.052 ± 0.002 


2.07 ±0.11 


-0.450 ±0.152 


20 


0.045 ± 0.002 


2.22 ± 0.09 


-0.556 ±0.133 


22 


0.040 ± 0.001 


2.36 ± 0.08 


-0.689 ±0.140 


24 


0.033 ± 0.001 


2.34 ± 0.06 


-0.681 ±0.127 



fitted very well by the function 

A N 



W (A) = 



In < A — e^v 



(A-e N ) 2 -l 



1/2 



B N , 



(18) 

where ejv, An and Bn are fitting parameters. As shown 
in Fig. [BJ all the data for N = 16, 18, 22, 24 and A < 5 
fit very well to Eq. (fl~8"|) . The results of these fits are 
collected in Table HIl 

To analyze the finite-size dependence in more detail, we 
adopt the standard finite-size scaling hypothesis^. We 
assume that in the infinite system, the DW width plays 
the role of the correlation length, that is, we assume that 



W (A) - W Q (A - 1)- 



(19) 



where v is a critical exponent. Finite-size scaling predicts 
that the effective critical value A^ = 1 + ejv where €n 
is proportional to N^ 1 ^ . Taking v = 1/2, Fig. [7] shows 
that A N converges to one as N increases. 



IV. THE STABILITY OF DOMAIN WALLS 

To describe the stability of the DW structure, we in- 
troduce <5„ (A) (n = 1, 2, N): 



Sn (A) = vWi^A)] 2 -^,^) 2 , (21) 



where 



_ £ {S* n {t)f dt 



[S'fa.tajA)]' 



(22) 



In order to show the physical meaning of 5 n , we write 

(SnW)as 



(s z n (t)> = c n + n n (t) , 



(23) 



where C n is a constant and fl n (t) is a time-dependent 
term. Then Eq. (|2"Tj) becomes 



t 2 -tl 



J t >n (t)dt 



2 \ 1/2 



(24) 



It is clear that if (S* (t)) is a constant in the time 
interval [ti,ta]» then 5 n (A) = 0. In general, since the 
initial state is not an eigenstate of the Hamiltonian Eq. 
([T]), the magnetization of each spin will fluctuate and 
^« (t) 7^ 0. If, after long time, the system relaxes to a 
stationary state that contains a DW, the magnetization 
of each spin will fluctuate around its stationary value 
C n . The fluctuations are given by O n (t). If |0„ (t)\ is 
large, the difference between the actual magnetization 
profile at time t and the stationary profile C n may be 
large. From Eq. (f2~4")) . it is clear that 5 n (A) is a measure 
of the deviation of (S* (t)} from its stationary value C n , 
averaged over the time interval [ti, tz]- Thus, 8 n (A) gives 
direct information about the dynamics stability of the 
DW. 

Figure [TU] shows the distribution S n (A) for different 
values of A. We only show some typical results, as in Fig. 



FIG. 6: The DW width as a function of A in a spin chain of N = 16, 18, 22, and 24 spins. The black dots are the simulation 
data and the solid line in each panel is given by Eq. (|18p . 



[3J As expected, the distribution of 5 n (A) is symmetric 
about the centre of the spin chain (n = 10.5). 

We first consider how S n (A) changes with A for fixed 
n. From FigfTUl we conclude: 

1) For the spins which are not located at the DW cen- 
tre, i.e., n ^ 10, 11, S n (A) decreases if A becomes larger. 
This means that the quantum fluctuations of these spins 
become smaller if we increase the value of A. This is 
reasonable because with increasing A, the initial state 
approaches an cigenstate of the Hamiltonian for which 
S n (A) =0 (Ising limit). 

2) For the spins at the DW centre, i.e., n = 10,11, 
when A becomes larger and larger, 6 n (A) first increases 
and then decreases. Qualitatively, this can be understood 
in the following way. When A is close to 1, the magneti- 
zation at the DW centre disappears very fast and remains 
zero. However, if A >> 1, the magnetization at the DW 
centre will retain its initial direction, hence the behavior 
of the spin at the DW centre will qualitatively change 
as A moves away from the critical point A = 1. In Fig. 
[TTl we plot Sio (A) (= <5n (A)) as a function of A. It is 



clear that <5io (A) first increases as A increases, reaches 
its maximum at A — 1.3, and then decreases as A be- 
comes larger. 

Now we consider the n-dependence of 5 n (A) for fixed 
A. Since S n (A) is a symmetric function of n, we may 
consider only one side of the whole chain, e.g., the spins 
with n — 1,2,..., N/2. From FigrTUJ according to the 
value of A, there are three different regions: 

1) 1.05 < A < 1.3: starting from the boundary 
(ri = 1), S n (A) first decreases, then increases, and finally 
decreases again as n approaches the DW centre (n = 10). 
As we discussed already, the fluctuation of the magneti- 
zation at the DW centre is small when A is close to 1. 
The spin at the boundary only interacts with one nearest 
spin, so it has more freedom to fluctuate. For the others, 
because of the influence of the DW structure (or bound- 
ary), the fluctuations of the spins which are near the DW 
(or near the boundary) are larger compared to those of 
a spin located in the middle of a polarized region. Thus 
S n (A) is larger if the spin is located near the DW or near 
a boundary. 
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FIG. 7: Fit of A* N to A* + A ■ N with A* = 1.009 ± 0.002, 
and A = 14.253 ± 0.660. 




FIG. 8: The DW width as a function of A (1.06 < A < 20) 
in the ground state of subspace M = in a spin chain of 
N = 20 spins. The black dots are the simulation data and 
the solid line is given by Eq. JTSJl, with ejv = 0.010 ± 0.001, 
A N = 1.87 ± 0.04 and B N = -0.550 ± 0.079. 



FIG. 9: The exponent C as a function of A max in a spin chain 
of N = 20 spins. The exponent C obtained by fitting the DW 
width to Eq. (|20|l . with A^ =2 o as obtained from the fit shown 
in Fig. [7] for A in the range [1.06, A max ]. 
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2) 1.3 < A < 5: & n (A) reaches its maximum at the 
DW centre. The reason for this is that in this regime the 
magnetizations of all spins retain their initial direction, 
therefore the spins that are far from the centre fluctuate 
little. 

3) 5 < A: in this regime (Ising limit), the initial state 
is very close to the eigenstate, and the fluctuations are 
small, even for the spins at the DW. 



V. SUMMARY 

In the presence of Ising-like anisotropy, DWs in a fer- 
romagnetic spin 1/2 chain are dynamically stable over 



FIG. 10: (Color online) 5 n (A) as a function of n for different 
A. Here ii = 101t, = 200r. We only show the data for 
A = 1.05, 1.06, 1.1, 1.2, 1.3, 1.5, 2, 5 and 20. The total number 
of spins in the spin chain is N = 20. 



extended periods of time. The profiles of the magnetiza- 
tion of the DW are different from the profile in the ground 
state in the subspace of total magnetization M = 0. 
As the system becomes more isotropic, approaching the 
quantum critical point, the width of the DW increases as 
a power law, with an exponent equal to 1/2. 
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FIG. 11: 5io (A) as a function of A. Here ti = lOlr, ti = 
200r. The total number of spins in the spin chain is N = 20. 
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